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ABSTRACT 

The  problem  of  finding  the  laminar,  compressible  boundary  layer  flow 
over  a  rotating  cone  at  small  angles  of  attack  is  considered.  The  cone  angle 
and  free  stream  Mach  number  are  such  that  the  non- viscous  flow  is  conical. 

Also  a  Prandtl  number  unity,  and  a  linear  viscosity- temperature  relation  are 
assumed.  The  "exact"  solution  for  small  angle  of  attack  and  small  spin  is 
obtained  by  a  perturbation  method  including  the  first  order  interaction  effects . 
For  a  slender  cone  the  results  are  used  to  calculate  the  Magnus  effects  due  to 
displacement  thickness.  As  yet  there  are  no  experimental,  measurements  which 
yield  a  conclusive  test  of  the  analytical  results  for  Magnus  force  given  here. 
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=  coordinate  angle  around  circumference  of  cone 
=  angular  velocity  of  cone  about  axis 
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=  indicates  reference  condition  taken  a&  conditions  on  surface 
of  cone  for  zero  angle  of  attack. 

=  subscript  indicating  free  stream  value. 

=  subscript  indicating  outer  edge  of  boundary  layer.. 

=  indicates  quantity  with  dimensions,  except  in  Eqs.  (l) 
and  (6)  -  (19). 
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INTRODUCTION 


The  study  of  three-dimensional  boundary  layer  flows  is  becoming 
increasingly  important.  The  class  of  flows  for  which  "exact"  solutions 
can  be  obtained  is  limited  to  flat  plate  and  yawed  cylinder  type  problems.^ 
However,  if  one  considers  flows  for  which  the  departure  from  a  basic  two- 
dimensional  flow  is  small,  then  by  application  of  small  perturbation 

1  2 

methods  a  much  wider  class  of  flows  becomes  .amenable  to  calculation.  * 

It  is  interesting  to  note  that  the  boundary  layer  flow  over  a  rotating  ' 
body  of  revolution  at  zero  angle  of  attack  can  be  classed  roughly  between 
the  two-and- three-dimensional  problems.  For  small  relative  spin,  however, 
this  problem  is  fairly  easy  to  solve. ^  A  body  of  revolution  with  non¬ 
zero  angle  of  attack  which  is  also  rotating  is  a  rather  intricate  three- 
dimensional  problem  but  one  which  is  of  basic  importance  in  ballistics. 

To  the  author's  knowledge  only  one  example  of  this  type  of  flow  has 
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been  considered.  Martin  has  solved  the  case  of  a  semi-infinite  open- 
ended  cylinder  in  incompressible  flow.  The  main  purpose  of  his  work  was 
to  test  a  hypothesis  that  Magnus  effects  on  slender  bodies  at  small  angles 
of  attack  are  caused  by  the  displacement  effect  of  the  boundary  layer. 

Since  the  prediction  of  Magnus  effects  has  remained  rather  elusive,  the 
fact  that  Martin  obtained  some  meaningful  results  was  encouraging. 

However,  because  of  the  model  chosen,  experimental  verification  could 
not  be  conclusive.  Thus  it  was  decided  to  try  a  similar  analysis  for 
a  more  easily  realizable  model.  Since  it  is  desirable  to  take  into 
account  compressibility  for  ballistic  applications,  the  cone  in  super¬ 
sonic  flow  seems  to  be  the  obvious  choice  for  an  analytical  aolution 
including  angle  of  attack  and  spin  effects.  This  is,  of  course,  because 
when  the  external  flow  is  conical  there  is  no  pressure  gradient  along 
generators,  and  this  allows  a  boundary  layer  solution  to  be  obtained  by 
introducing  similarity  variables. 

The.  three  special  cases  of  the  above  problem  have  been  worked  out; 
for  zero  angle  of  attack  and  zero  spin  by  Hantsche  and  Wendt p  for  zero 
spin  by  Moore and  for  zero  angle  of  attack  by  Illingworth.^  Moore's 
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analysis  is  restricted  to  the  first  power  in  angle  of  attack,  and  Illingworth 
considers  up  to  cubic  terms  in  the  spin.  If  one  is  interested  in  only  first 
order  effects  of  both  angle  of  attack  and  spin,  these  two  results  can  be 
superposed  since  the  perturbations  enter  linearly,  However,  to  obtain  the 
interaction  effects,  which  are  necessary  for  Magnus  effects,  additional 
analysis  must  be  made.  This  is  what  has  been  done  here.  Only  laminar  flow 
is  considered;  compressibility  is,  of  course,  not  neglected;  heat  conduction 
is  allowed  for  with  the  assumption  of  Prandtl  number  unity  and  a  linear 
viscosity- temperature  relationship.  In  the  solution  to  the  boundary  layer 
problem  the  only  other  assumption  (aside  from  the  representation  in  terms  of 
small  perturbation)  is  that  the  external  flow  Mach  number  and  the  cone  angle 
are  such  that  the  flow  is  conical.  When  the  results  are  applied  to  Magnus 
force  calculations  slender  body  theory  is  used  so  that  the  restrictions  to 
small  cone  angles  and  slightly  supersonic  Mach  numbers  are  necessary. 

Aside  from,  the  interest  in  Magnus  effects  the  three-dimensional  problem 
considered  here  has  some  intrinsic  interest  so  that  additional  results  are 
given,  e.g.,  velocity  profiles  and  skin  friction.  In  particular,  the 
following  should  be  noted.  To  obtain  generally  useful  results  for  three- 
dimensional  flows  approximate  (momentum-integral)  methods  must  he  developed; 
this  is  especially  true  if  the  turbulent  case  is  to  be  considered.  Reliability 
of  these  methods  cannot  be  based  alone  on  experience  in  two-dimensional  flows 
since  there  are  fundamental  differences  in  the  two  cases.  Thus  exact  solutions 
are  necessary  to  test  proposed  approximate  methods. 

The  Magnus  effects  to  be  considered  here  are  those  commonly  observed  in 
ballistics.  In  the  terminology  of  airplane  dynamics  the  Magnus  effects  are 
a  side  force  and  yawing  moments  due  to  angle  of  attack  which  implies  that 
the  body  must  be  spinning  for  non-zero  Magnus  effects.  It  is  clear  that  non- 
viscous  flow  could  not  allow  Magnus  effects  (this  is  discussed  more  fully  in 
Reference  4)  thus  one  is  forced  to  investigate  the  boundary  layer  flow. 
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Assuming  a  Prandtl  number  equal  to  one  and  a  linear  viscosity- temperature 

relation  give  the  boundary  layer  solution  for  no  heat  transfer  at  the  surface 

only  if  the  surface  is  stationary.  For  a  spinning  body  this  is  not  the  case, 
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as  Illingworth  has  noted.  For- the  cone  one  finds  that,  exactly,  the  solu¬ 
tion  is  for  a  surface  temperature  which  varies  as  the  square  of  the  distance 
along  a  generator.  An  alternate  interpretation  is  that  the  heat  transfer  at 
the  surface  is  zero  neglecting  quadratic  terms  in  a  and  (Mr.  Martin  Fiebig 
called  to  the  author’s  attention  an  error  in  the  original  manuscript  concern¬ 
ing  this  point.) 

EQUATIONS  OF  MOTION 

The  form  of  the  equations  of  motion  used  here  is  taken  from  the  work  by 
8  *7 

Moore  *  and  will  only  be  summarized.  A  two-component  vector  potential  is  in¬ 
troduced  to  satisfy  the  continuity  equation  exactly  (generalization  of  stream 
function) .  Prandtl  number  unity  and  the  viscosity  temperature  relation 

H  /  h  =  CT  /  T 

are  assumed,  where  the  bar  refers  to  the  reference  condition  taken  as  the  non- 
viscous  flow  at  zero  angle  of  attack  evaluated  at  the  surface  and  C  is  a 
suitably  chosen  constant.  Two  transformations  are  introduced;  one  like  that 
of  Howarth  to  try  to  remove  compressibility  effects  and  another  like  that  of 
Mangier  to  try  to  remove  curvature  effects.  Non- dimensional  dependent  variables 
are. used  with  the  reference  condition  denoted  by  a  bar  with  the  same  meaning 
as  above,  and  lengths  are  made  dimensionless  with  the  length  L  =  (Cfi)  /  (p  u). 
Let  u  and  w  denote  (non-dimensional)  velocities  along  generators  (x  coordinate) 
and  around  the  circumference  (0  -  coordinate);  see  Figure  1.  If  F  and  G  denote 
the  two  components  of  the  vector  potential  mentioned  above  then 

u  =  F  ,  w  =  G 

v  n 

where  the  subscript  denotes  differentiation  and 

n  =  x(p  /  p)1/2  f  pdy 
o 
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Let  0  be  the  half- angle  of  the  cone 


0  =  sin  0  j  |  =  x  /3 

then  the  equations  of  motion  are  (Equation  Reference  7). 

'  [V  V{J8,3  Fnn  +  0iV(585)  ’  (G//(3s) 

P’(0)  GF  / (60gp)  +  F 


w 


Vie  ■  ft  +  V(59^  am  *  <VV(5e5)  *  W<*>  - 

*-  P’  (0)/(5©6p)  +  P'(0)GG^/(69|p)  +  GTJTlT1 


T  +  (F  )  +  (G  )2  =  Tx  +  u^  + 


(la) 


(lb) 


(lc) 


where  p,  p,  and  T  are  (non-dimensional)  pressure,  density,  and  temperature j 
p'(0)  is  the  derivative  of  pj  and  the  subscript  1  denotes  conditions  at 
the  outer  edge  of  the  boundary  layer.  The  form  of  the  variables  7  and  | 
comes  from  the  Hovarth  and  Mangier  transformations.  For  the  boundary 
conditions  at  the  outer  edge  of  the  boundary  layer  (y  =  tj  =  CD  ). 


u  =  u^  =  1  -  aA^  cos  0 


w  =  wn  =  sin  0 

1  (2) 

p/p  =  1  +  aAj  cos  0 

P1  =  1  +  aA^  cos  0 
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since  the  flow  is  conical.  The  A^  are  tabulated  ,  see  Reference  6  for 
further  discussion.  It  is  assumed  that  the  rotation  of  the  cone  does  not 
alter  the  external  flow  so  that  the  pressure  is  given  by  the  expression 
in  (2).  For  the  boundary  conditions  on  the  body  surface  (y  =  t}  =  0) 
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u  »  v  ■  0 


V  ■  (BT  /  U1  (3) 

r  ■  xL6 

where  m  is  the  angular  velocity  of  the  cone  about  its  axis,  r  is  the 
local  radius  and  L  is  reference  length  defined  above. 

Since  the  pressure  gradient  is  zero  along  the  generators  of  the  cone 
a  Blasius  type  similarity  variable  is  introduced. 

».  - 1/(2  vD- 

(This  differs  by  a  factor  of  one -half  from  the  X  used  by  Moore.)  However, 
because  of  the  spin  of  the  body  there  is  no  longer  similarity  in  terms 
of  X.  Thus  variation  with  x  must  be  allowed,  and  for  this  it  is  convenient 
to  introduce  a  new  variable  *  , 

<  =  £DT  /  u~  =  cdL6  (3|  /  u* 

Furthermore  the  following  forms  for  F  and  G  are  assumed: 

F  =  |1/2  f ( x  ,  X,  0) 

G  =  |1/2  g(x  ,  X,  0) 

v 

For  most  ballistic  applications  the  maximum  vqlue  of  fc  will  be  less  than 
0*5>  but  it  is  assumed  here  that  x  remains  small  enough  so  that  f  and  g 
can  be  expanded  in  power  series  in  x  .  Hence,  for  this  analysis,  both 
the  length  of  the  cone  and  the  angular  velocity  are  limited.  Thus  we 
write 

f  =  fQ(x,  0)  +  x  f-jU,  0)  +  ...  . 

g  =  gQ(x,  0)  +  x  g^X,  0)  +  .  .  .  (h) 

T  =  tQ(x,  0)  +  x  tx(x,  0)  +  .  .  . 
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For  u>  =  0,  fg  and.  g^  are  the  same  as  the  f  and  g  defined  by  Moore  .  For 

a  =  0  all  quantities  are  independent  of  0,  and  we  obtain  essentially 

3 

Illingworth's  first  order  spin  effect  terms ;  the  correspondence  is 

(Vj  -  ro 

(eO  -  s,  /2  . 

I  X 

where  the  I  indicates  Illingworth's  functions  (f^  would  be  identically 
zero  for  a  =  0). 

The  partial  differential  equations  satisfied  by  the  coefficient  of  k 
in  (4)  will  not  be  written  down.  In  order  to  reduce  the  integration 
problem  to  one  amenable  to  numerical  processes  the  dependence  of  the 
coefficients  on  0  is  assumed  as  follows  (the  form  can  be  deduced  from  the 
appropriate  boundary  conditions): 

fQ  =  FqqCx)  -  a/^  cos  0  fq1U)  +  .  .  . 

gQ  -  0^2  sin  0  VX>  +  *  ‘  ’ 

(5) 

f,  =  sin  0  F^Cx)  +  .  .  . 

&L  =  Gl0(^)  "  cos  0  outo  +  •  ■  • 

Only  the  functions  of  X  explicitly  given  in  the  expansions  (5)  will  be 

obtained.  If  the  assumed  form  of  the  solution  be  regarded  as  a  power 

series  expansion  in  the  two  small  parameters  a  and  k  it  is  evident  that 

only  one  second  order  term  is  included,  i.e.,  the  product  (interaction) 

2  2 

term  a  k  .  The  a  and  k  terms  could  be  included  at  the  expense  of  more 
numerical  integrations.  This  may  seem  inconsistent,  but,  at  least  for 
Magnus  effects,  the  second  order  terms  are  unimportant .  The  functions  of 
X  indicated  in  (5)  can  be  identified  as  follows;  F^  governs  the  pure 
axial  flow,  FQ1  and  GQ1  the  angle  of  attack  effects,  G1Q  the  spin  effect, 
and  F11  and  G_L1  the  interaction  effects. 
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As  expected  F_q  satisfies  the  Blasius  equation.  The  remaining 
functions  satisfy  linear  third-order  ordinary  differential  equations 
with  boundary  conditions  at  X  =  0  and  X  =  CO  •  These  were  not  the 
equations  that  were  integrated,  however,  and  therefore  will  not  he  written 
down,  for  the  following  reason.  These  differential  equations  contain, 
through  the  A^,  the  "Wo  parameters  cone  angle  ©  and  free  stream  Mach 
number  Mq.  But  since  they  are  linear  it  is  an  easy  matter  to  rewrite  the 
equations  in  a  form  which  contains  universal  functions,  i.e.,  independent 
of  any  parameters.  Thus  the  following  sequence  of  equations  was  obtained 
(the  Blasius  equation  and  that  for  G^q,  which  is  related  to  Illingworth's 
g^,  require  no  change,  but  are  rewritten  in  a  different  notation): 

Y"  +  •  =  0 

Vo)  =  V(0)  -  0,  Y  (  CO  )  =  2  (6) 


F00  ■  *1 

Y,f  +  Vi"  -  (2/3) YV  =  -  (fi/3)d  -  \'2/k) 

Vo)  -  Y  (0)  =  o,  Y  (©  )  =  o  (7) 

Then  G0l  =  hj  +  q  q  =  1  +  (7  -  1)mP  /2 


Define : 

Vf)  =  f  * 1  +  y’  1  +- Y  ,f 

VV  =  0  h^O)  =  ^'(0)  =  0  (8) 

b2'  (  00  )  =  1 

VV  -  h1,,k1  =  0  h?(0)  =  h3«(0)  =  0  (9) 

h  '  (  OO  )  -  0 

3 


■x 
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Then  r  .  1 

F01  =  [2^/(5©^)]  h-L  +2[l  -  (2^/59^)]  Hg  +  [2Ag  q/i&Aj]  h? 

Define : 

Lg(f )  =  f  •  • '  +  l^f  "  -  (4/5)  f  • 

W  -  0  kg(0)  =  0,  kg* (0)  =  2, 

kg'  (oo  )  =  0 

(10) 

Then 

G10  =  k2 

/ 

Also: 

I^)  +  kg"  hg  -  (4/^)k2l  hg*  =  0 

(H) 

L2(k^)  +  k2»»  (h3  -  kx)  +  k2'  [kx«  - 

(4/3)  h5’  ]  =0 

(12) 

I^)  -  (1/5)1^  k2«  =  0 

(15) 

\(0)  =  k± ' (0)  =  k1«  (GO  )  =  0, 

i  -  5,4,3 

Then 

G11  =  [2/(^0)]  {[2(30A1/  2Ag)  -  l] 

k3  +  qk^  +  k5| 

Define : 

L  (f)  =  f"  +  hxf "  ■  -  (2/3)h1*  f  + 

(5/3)  Y*  f 

L5(h4)  +  (4/3)  kg'  h1'  =  0 

(14) 

L^)  +  (4/3)  kg'  k”  =  0 

(15) 

.  L5(h6)  -  hx'  kg'  +  h1"  k5  =  0 

(16)  . 

L3(V  '  V  V  +  V  k4  =  0 

(17) 

W  ’  V  hg'  +V'k5  =  0 

(18) 

L3(V  '  k2  hl,,=0 

h^o)  =  lyto)  =  ht‘  ( co  )  -  o. 

i  =  4,5,  .  .  .,  9. 

(19) 

15 


Then 


Fn  -  (Ag/V^  +  qV  +  (2/36)2  (A2/AL)(h6  +  qV  +  CV30)(1  -2A2/36A1)h8 
+  (Ay^)^ 

Equations  (6)  -  (19)  were  integrated  numerically  on  the  QRDVAC  and  are 

tabulated  in  the  Appendix.  The  functions  were  obtained  to  an  accuracy 

of  at  least  five  significant  figures.  More  accuracy  was  needed  .in  the 

first  few  equations  because  of  accumulation  of  errors.  Equations  (7) , 

6 

(8)  and  (9)  are  essentially  the  same  as  those  integrated  by  Mooreu  (the 
numerical  results  agreed  except  that  for  the  larger  values  of  X  the 
agreement  was  only  in  the  second  or  third  significant  figure).  Equation 
(10)  is  essentially,  the  same  as  one  of  Illingworth’s  (results  agree  to  as 
many  places  as  given).  The  remaining  equations  give  the  interaction 
effects . 


RESULTS 

After  numerical  integration  of  the  sequence  of  ordinary  differential 
equations  it  is  possible  to  compute  a  number  of  quantities  of  interest. 

(a)  Velqcity  Profiles: 

The  non-dimensional  velocity  components  u  and  w  are  obtained  from 

-  V  V2 
•  v  •  S  =V2 

and  making  use  of  (4)  and  (5).  For  given  free  stream  Mach  number,  M, 
and  cone  angle,  0,  ,  the  velocities  are  functions  of  a ,  k,  and  0. 

In  Fig.  2  u  and  v  are  plotted  for  the  case  M  =  1.82,  ©  =  10°,  a  =  1°, 

K  =  0.1  as  functions  of  X  for  0  =  -  n/2,  0,  it/2,  and  it. 
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FIGURE  2 

Non-dimensional  velocity  components  u,  w  for  M  *  1.82,  0=10°,  a  *  1°,  *  *  0.1 


As  expected  there  is  a  reversal  in  the  w  component  for  0<O,  i.e. , 
where  the  rotational  velocity  opposes  the  direction  of  the  flow  outside 
the  boundary  layer.  For  the  values  of  parameters  chosen  the  u  component 
does  not  change  drastically  around  the  circumference.  There  is  no 
indication  of  separation  type  profiles.  Note  that  it  has  not  been  necessary 
to  compute  the  temperature  variations.  However  this  would  have  to  be. 
done  if  the  velocity  profiles  were  derived  as  functions  of  the  physical 
coordinate  y. 


(b)  Skin  Friction: 

The  two  components  of  skin  friction,  expressed  in  dimensionless  form 

are 


where  the  prime  indicates  that  a  quantity  has  dimensions.  To  the  same 
order  of  approximation  that  has  been  carried  throughout,  the  following 
expressions  are  obtained. 


2cfx  =  C*/*)172  Foo"(0)  ■  a{AlF-Ol’,(0)  “  U/2)A^00,,(0)}cos  0 


+  K 


a  A^F^iO)  sin  0 


(20) 


2Cf0  =  W*x)1/2  [*Glo"(0)  +  a  A2G0i',(0)  Sln  * 

+  *  a  |a5G10m  (0)/2  -  A2G11m(0)  ■  cos  0] 


(21) 


where  Rx  =  x'  p'  u'/  p* .  The  terms  in  (20)  and  (21)  which  are  independent 
of  or  linear  in  a  and  K  can  be  obtained  from 'References  6  and  3  respectively. 
The  interaction  terms  (a  k  )  ore 
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a%*(V*>1/2/^  a  sin  0 =  ^11"  (°)/2 

=  [.6145  +  .2522  T  -1  +  (2/56)2  (.3592 

-  (2/36)  [.5109  \  +  .2395 

.aVV5 c)1/2/s  *  s  a  cos  0  -  E»  V<°)  '  2A20ll“(°>]  !k 
=  .4896  A^_  -  .4898  -  (A2/36)  •  [..4496  +  .5068 

where  T  ^  =  (7  -  l)M^/2.  The  right  hand  sides  of  equations  (22)  and  (23) 
are  plotted  in  Figures  3  and  4  respectively  for  cone  angles  of  10°,  12. 5°, 
and  15°^  According  to  the  expression  (20)  for  the  effect  'Of  the  a  k 
term,  for  negative  0,  is  to  reduce  the  skin  friction.  This  reduction 
is  a  maximum  at  0  =  -«/2.  However,  for  resonably  small  values  of  a  and 
k  the  reduction  is  not  enough  to  indicate  component  separation.  For 
example  to  reduce  Cf  to  zero  for  =  10°  at  0  =  -  it/2  it  would  be 
necessary  to  have  a  k  -  .06.  For  such  values  of  a  and  k  the  higher  order 
terms  would  no  longer  be  negligible.  Since  this  is  a  three-dimensional 
flow  it  is  by  no  means  clear  that  component  separation  is  pertinent  to 
the  general  question  of  separation. 

—  l/2  ”  *  ♦ 

In  Fig.  5  &  sketch  of  C„j(R  /3C)  is  presented  for  M  =  1.82  and 
~  o  1)0  x 

tv  =  10  .  It  is  seen  that  the  variation  of  the  skin  friction  is  affected 
significantly  by  the  interaction  terms.  However,  the  integrated  value, 
which  gives  the  spin  reducing  torque,  is  unaffected  in  the  present 
approximation .  (See  Reference  3  for  a  discussion  of  the  effect  of  the 
second  order  term  in  the  spin  on  this  torque . ) 

(c)  Displacement  Thickness 

To  find  the  displacement  thickness  for  a  three-dimensional  boundary 
layer  flow  it  is  necessary  to  solve  a  first  order  partial  differential 
equation"^1 For  the  cone  problem  if  A  is  the  displacement  thickness 
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3gCfx(Rx/3C) 
dKda  sin  4* 


12 
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dKdCL  cosrf) 


2: 


FIGURE  5 


VARIATION  OF  Cfj(I ^/jC)1/2  AROUND  CIRCUMFERENCE  FOR  M  -  1.82  andOfr-  10°. 


L  ■  C*  jl  /  p  u 


Integrating  (24)  and  introducing  a  displacement  thickness  Reynolds  number 
ra  =  p*  uvy  p* 
the  following  is  obtained 


where  1^,  and  are  functions  of  free  stream  Mach  number  and  cone 
angle.  From  References  10  and  6,  and  can  be  obtained  and 


J2  =  Kk  +  +  (Kg/5©) 

where  =  (2/50)  jA^  (.41 42  -  .5198  T  _1)  -A^  (.2158  +  .2771  T 
+  A2  £.4106  +  .6396  T  _1  +  .1697  T  "2 
+  (2/3© f  (.3186  +  .5247  T  "1  +  .1495  T  "2)1 


=  A2  T  “1  (.6692  +  .2495  I  _1) 

K6  =  a2  (2/30) (.9207  +  3289  T  “1)  -  .6546  -  .6558 


The  a  *  term  in  (25)  makes  the  displacement  thickness  unsymmetrical  with 
respect  to  the  plane  of  yaw  and  is  indicative  of  Magnus  effects.  For  these 
effects,  which  will  be  considered  in  the  next  section,  it  is  convenient 
to  have  J2U//u  which  is  plotted  in  Fig.  6  where  U  is  the  free  stream 
velocity.  This  velocity  ratio  is  introduced  because  the  dimensionless  force 
and  moment  coefficients  are  conventionally  based  On  free  stream  velocity. 
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u/n 


M 

FIGURE  6 

CONTRIBUTION  OF  INTERACTION  TERM  TO  DISPLACEMENT  THICKNESS  A. 
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In  Fig.  T  the  variation  of  A  around  the  circumference  is  shown  for  M  =  -1.82, 

10°.  The  effect  of  the  spin  on  displacement  thickness  is  to  rotate 
it  in  the  direction  of  the  spin.  For  the  value  of  a  and  k  of  Fig.  7  the 
maximum  and  minimum  are  rotated  through  15°  from  the  plane  of  yaw. 

* 

MAGNUS  EFFECTS 

The  concept  that  the  unsymmetrical  displacement  thickness  on  a 
rotating  body  at  angle  of  attack  causes  Magnus  effects  was  used  by 

4 

Martin  for  the  case  of  a  semi-infinite  cylinder.  Basically  the  same 
idea  will  be  used  here.  The  potential  flow  over  the  distorted  body,  i.e.r 
body  plus  displacement  thickness,  is  calculated  by  means  of  slender  body 
theory.  Thus  we  now  assume  that  the  cone  angle  is  small.  The  displacement 
thickness  as  given  by  (25)  is  measured  normal  to  the  body  st^rface  but  for 
the  approximation  of  slender  body  theory  this  thickness  can  be  added  to 
the  cross-section  of  the  body  normal  to  its  axis.  Then  the  cross-section 
of  the  bbdy  can  be  expressed  in  polar  coordinates  as 

r'/Z  =  a  +  b  cos  0  +  c  sin  0  (;26) 

where  a  =  (x1  ®:/&)  +  2K1(Cx« 
b  =  2a  J^Cx’/SR/)1'2 
c  =  -2a  J2£l(x'/i)(Cx'/5R/)1//2 

ft  =  CD  @  jj/  U* 

with  l  the  length  of  the  cone,  the  Reynolds  number  based  on  l  and  the 
surface  quantities  at  zero  angle  of  attack,  and  the  prime  indicates  a 
quantity  with  dimensions.  To  the  first  order  in  a  -the  expression  (26) 
for  the  cross-section  is  a  circle  of  radius  a  and  center  (b,  c). 

The  details  of  the  slender  body  calculation  need  not  be  given  since 
12 

Ward  has  presented  all  the  necessary  formulas  for  the  cross  force  and 
moments.  Ward  uses,  instead  of  the  axis  of  the  body,  the  s-axis  parallel 
to  the  direction  of  the  free  stream.  With  respect  to  the  s  -axis  the 
center  of  the  circle  (26)  is  (b*,  c)  Where 
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FIGURE  7 


VARIATION  OF  (R^/2C)(x/3)^^  AROUND  CIRCUMFERENCE  FOR  M  -  1.82  and  (Si  =  10 


b*  =  b  +  as 


and  for  the  approximation  made  here  x:  =  s .  Then  the  complex  transverse 
force,  Fj  is  given  by 

.  *  v 


2F/p  U^S  =  -  2(db*/da  +  idc/ds)  _  ,  (2j) 

O  S  —  X 

* 

where  S  is  the  base  area  and  pQ  the  free  stream  density.  The  real  and 
imaginary  parts  give  the  lift  and  cross-force  respectively.  The  negative 
sign  in  (27)  results  from  the  choice  of  axes  in  the  cross-section.  Thus 
the  displacement  thickness  caused  a  decrease  in  the  lift  coefficient  from 
the  value  2a  by  an  amount 


2aJ1(c/^l)1/2 


(j^  is  negative).  The  Magnus  force  coefficient  (in  aerodynamic  notation) 
is 


=  6  OJ^C/SF^) 


1/2 


(28) 


The  ratio  of  the  Magnus  force-  to  lift  force  due  to  displacement  thickness 
is 


%  J2/J-L 

which  for  ©  =  10°  is  of  the  order  of  magnitude  10ft  which  is  less  than  one 
for  the  small  spin  contemplated  here.  Such  a  small  Magnus  force  would  be 
difficult  to  measure. 

The  center  of  pressure  of  the  Magnus  force  is  found  to  be 

1  *f  8 ^  [f  +  2VCS/5R//2]  ds/  r@;  +  2K1(C/3Ri)1/^j 

in  units  of  the  length  of  the  cone.  To  within  a  few  percent  at  most 
this  result  shows  that  the  center  of  pressure  is  independent  of  Mach 
number  and  Reynolds  number  and  is  given  by 

C.P.  =  (5/7) i. 
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In  ballistic  notation  the  Magnus  force  and  moment  coefficients,  Kp 
and  respectively,  are  found  to  be’ 

Kp  =  F/  p^d2  V  a  =  («/Q)(J2v/u')(xM£)l/2 

^  =  T/  p^d3  V  a  =  -(C.P.J,  -  C.MOKp 

where  d  is  the  diameter  of  the  base,  v  =  cod/u,  C.P.  is  the  center  of 

X 

pressure  and  C.M.  is  the  center  of  mass  both  measured  in  calibers 
(diameters).  The  moment  is  taken  about  the  center  of  mass,.  In  Pig.  6 
J^U/ u'  is  plotted. 

As  yet  there  appear  to  be  no  experimental  measurements  which  provide 

a  conclusive  test  of  these  analytical  results.  There  is  some  work  in 

progress  at  the  BRL  which  may  yield  such  a  test.  Some  measurements  have 
13 

been  reported  giving  only  the  Magnus  moment  on  a  slender  cone.  However, 
the  boundary  layer  was  definitely  turbulent  over  most  of  the  cone. 
Moreover,  it  is  felt  that  the  predicted  moment  will  be  more  in  error  than 
the  force.  When  the  above  result  for  K^,  was  compared  with  the  results  of 
Ref.  13  it  was  found  to  disagree  considerably  but  in  view  of  the  above 
statements  this  is  not  regarded  as  conclusive. 


Finally  the  contribution  to  the  Magnus  force  of  the  meridional  skin 

friction  will  be  discussed.  This  was  done  for  Martin's  cylinder  problem 
14 

by  Kelly  who  found  that  skin  friction  contributed  a  negative  Magnus 
force  about  7  percent  of  displacement  thickness  contribution.  Integrating 
(21)  it  is  found  that  for  the  cone  the  a  k  term  gives  a  negative  Magnus 
force  which  for  ©  =10°  is  about  20  percent  of  that  due  to  displacement 
thickness. 

4 

Martin  has  attempted  to  extend  his  analysis  of  the  laminar  flow 
over  a  cylinder  to  the  turbulent  case  by  assuming  that  the  Magnus  force 
coefficient  depends  on  the  product  of  displacement  thickness  at  the  base 
for  zero  angle  of  attack  and  the  length  of  the  cylinder  as  in  the 
laminar  case.  Aside  from  the  fact  that  an  empirical  constant  is  left 
to  be  determined  it  is  not  clear  that  the  above  assumption  is  valid. 


* 
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Thus  the  Important  matter  of  treating  the  turbulent  boundary  layer  remains 
to  be  solved  as  well  as  the  extension  to  mare  general  bodies.  This  will 
necessitate  the  development  of  approximate  methods  for  three-dimensional 
boundary  layer  flows. 


f?.  S 

E.  SEENEY 
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Equation  7 
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1.001 67 

— 

.00756 

.9 

.39071 

.82808 

.71026 

3.4 

2.96469 

1.00105 

— 

.00500 

1.0 

.4  7693 

.89482 

.62291 

3.5 

3  .06478 

1.00065 

— 

.00322 

1.1 

.56937 

.95240 

.52766. 

3.6 

3.16483 

1.00039 

— 

.00203 

1.2 

.66708 

1.00021 

.42808 

3.7 

3.26486 

1 .00023 

."0124 

1.3 

.76908 

1.03801 

.32825 

3.8 

3.36488 

1.00013 

— 

.00075 

1.4 

.87436 

1.06599 

.23238 

3.9 

3.46489 

1.00007 

— 

.00044 

1  .5 

.98197 

1 .08475 

.14438 

4.0 

3.56489 

1.00004 

_ 

.00025 

1.6 

1 .09103 

1.09524 

.06750 

4.1 

3.66490 

1.00002 

— 

.00014 

1.7 

1.20078 

1.09870 

.00399 

4.2 

3.76490 

1 .00001 

.00007 

1.8 

1.31059 

1.09652 

- 

.04499 

4.3 

3.86490 

1.00000 

— 

.00004 

1.9 

1.41995 

1.09018 

- 

.0795? 

4.4 

3.96490 

1.00000 

— 

•000"2 

2.0 

1.52853 

1.08107 

- 

.10064 

4.5 

4.06490 

1 .00000 

— 

."0001 

2.1 

1.63611 

1.07044 

- 

.11019 

4.6 

4.16490 

1.00000 

— 

.00000 

2.2 

1.74260 

1.05934 

- 

.11047 

4.7 

4.26490 

l.ooooo 

— 

.000OO 

2.3 

1.84799 

1.04857 

- 

.10398 

4.8 

4.36490 

1 .00000 

— 

,00000 

2.4 

1.95235 

1.03869 

•09310 

4.9 

4.46490 

1*00000 

— 

*0000" 

Equation  9 

hi 

h1* 

\ 

hx 

5 

.0 

#nnno" 

#nrmoo 

— 

.38265 

2.5 

-  .42381 

-  .02969 

•09334 

.1 

-  .00191 

-  ."3825 

— 

.3821 ? 

2.6 

-  .42635 

-  .02145 

.07199 

.2 

-  .00764 

-  ."7632 

- 

.37850 

2.7 

-  .42816 

-  .01518 

.05413 

.3 

-  .01715 

-  .11375 

- 

.3691 0 

2.8 

-  .42944 

-  .01051 

.0397 n 

.4 

-  .03035 

-  .14987 

- 

.35169 

2.9 

-  .43031 

-  .00713 

.02843 

.5 

-  .04705 

-  .18377 

- 

.32472 

3.0 

-  .43090 

-  .00474 

.01987 

.6 

-  .06700 

-  .21447 

- 

.28739 

3.1 

-  .43128 

-  ,00308 

.01357 

.7 

-  .08981 

-  .24091 

- 

.23983 

3.2 

-  .43153 

-  .00196 

.00906 

.8 

-  .11501 

-  .26213 

- 

.18316 

3.3 

-  .43169 

-  .00123 

.00591 

.9 

-  ,  1  42.0  3 

-  .27731 

- 

.11945 

3.4 

-  .43179 

-  .00075 

,00377 

1  .0 

-  .17025 

-  .28588 

- 

.0516? 

3.5 

-  .43185 

-  ,0"045 

.00235 

1.1 

-  .19898 

-  .28761 

."1683 

3.6 

-  .43188 

-  .00026 

."0143 

1.2 

-  .22755 

-  .28262 

.08215 

3.7 

-  .43190 

-  .00015 

.00085 

1.3 

-  .25530 

-  .27140 

.1407? 

3.8 

-  .43191 

-  ,00"O8 

,0005" 

1.4 

-  .28165 

-  .25480 

.18948 

3.9 

-  .43192 

-  . 0" 004 

.00028 

1.5 

-  .30612 

-  .23391 

.22619 

4.0 

-  .4319? 

-  .oooo? 

.00016 

1.6 

-  .32833 

-  .21000 

.24972 

4.1 

-  .43192 

-  .00001 

,000"8 

1.7 

-  .34806 

-  .18441 

.26007 

4.2 

-  .43193 

-  .00000 

.00004 

1.8 

-  .36520 

-  .15839 

.25827 

4.3 

-  .43193 

-  .00000 

. 0  00"2 

1  .9 

-  .37976 

-  .13309 

.24621 

4.4 

-  .43193 

-  .0"00" 

,00001 

2.0 

-  .39187 

-  .10941 

.2263" 

4.5 

-  .43193 

-  .0"""" 

,000nr\ 

2.1 

-  .4017? 

-  .08801 

.20115 

4.6 

-  .43193 

-  ,0"OOO 

,0"00n 

2.2 

-  .40957 

-  .06927 

.17330 

4.7 

-  .43193 

-  .OOOOO 

2.3' 

-  .41567 

-  .05336 

.14497 

4.8 

-  .43193 

-  .000OO 

.Ooorm 

2.4 

-  *42033 

-  .04023 

♦  n  7  9 1 

4*9 

-  *43193 

-  *oonoo 

ionorin 
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Equation  10 

k"  ' 

\ 

.  k2 

*2 

.0 

.ooooo 

2. O0000 

-1.95932 

2.5. 

1  .30676 

.0041  0 

.1 

.  19021 

1.80464 

-1.94228 

2.6 

1.30710 

.00264 

.2 

.36103 

1.61257 

-1.89422 

2.7 

1.30731 

.00168 

.3 

.51293 

1.42668 

-1 .81951 

2  1 8 

1  .30744 

.00105 

.4 

.64666 

1.24941 

-1.72262 

2.9 

1  .30753 

.00064 

.5 

.76317 

1,08275 

-1.60808 

3.0 

1  .30758 

.00038 

.6 

.86361 

.92823 

-1 .48044 

3.1 

1.30761 

.00023 

.7 

.94926 

.78695 

-1 .34415 

3.2 

1  .30763 

.00013 

.8 

1.02147 

.65954 

-1.20350 

3.3 

1  .30764 

.00007 

.9 

1.08164 

.54626 

-1.06247 

3.4 

1.30764 

.00004 

1.0 

1.13118 

.44694 

-.  .92465 

3.5 

1  .30765 

.00002 

1.1 

1. 1 7148 

.36112 

.  -  .79308 

3.6 

1  .30765 

.00001 

1.2 

1.20383 

.28803 

-  .67024 

3.7 

1.30765 

.ooooo 

1.3 

1.22948 

.22672 

-  .55797 

3.8 

1  .30765 

.ooooo 

1.4 

1.24953 

.17605 

-  .45744 

3.9 

1.30765 

.ooooo 

1.5 

1.26500 

.13482 

-  .36924 

4.0 

1.30765 

.ooooo 

1.6 

1.27677 

.10179 

-  .29338 

4.1 

1.30765 

.00000 

1.7 

1.28559 

.07575 

-  .22940 

4.2 

1  .30765 

•  0n000 

1.8 

1.2  921 1 

.05554 

-  .17648 

4.3 

1  .30765 

.OOOOO 

1.9 

1.29686 

.04012 

-  .13356 

4.4 

1  .30765 

.OOOno 

2.0 

1.30026 

.02854 

-  .09941 

4.5 

1  .30765 

•ooooo 

2.1 

1.30267 

.01999 

-  .07275 

4.6 

1.30765 

.00000 

2.2 

1.30434 

.01378 

-  .05235 

4.7 

1  .30765 

-  .00000 

2.3 

1.30548 

.00935 

-  .03702 

4.8 

1.30765 

-  .ooooo 

2.4 

1.30625 

.00624 

-  .02573 

4.9 

1  .30765 

-  .ooooo 

Equation  11 


\ 

k’ 

5 

5 

5 

.0  • 

.ooooo 

.00000 

- 

.48957 

2.5 

— 

.32439 

— 

.01161 

.1 

- 

.00243 

- 

.04853 

- 

.47701 

2.6 

— 

.32537 

.00807 

.2 

- 

.00962 

- 

.09466 

- 

.44225 

2.7 

— 

.32604 

— 

.00550 

.3 

- 

.02122 

- 

.13639 

- 

.38972 

2.8 

— 

.32650 

— 

.00068 

.4 

- 

.03670 

- 

.17216 

- 

.32395 

2.9 

— 

.32680 

— 

.00242 

.5 

- 

.05542 

- 

.20089 

- 

.24950 

3.0 

— 

.32699 

— 

.00156 

.6 

- 

.0766? 

- 

.22192 

— 

.17089 

3.1 

— 

.32712 

— 

.00098 

.7 

- 

.09954 

- 

.23507 

- 

.09239 

3.2 

- 

.32  72  0 

— 

.00061 

.8 

- 

.12338 

- 

.24054 

- 

.01790 

3.3 

— 

.32725 

— 

.00037 

.9 

- 

.14741 

- 

.23890 

.04928 

3.4 

— 

.32727 

— 

.00022 

1.0 

- 

. 1 7095 

- 

.23101 

.10657 

3.5 

— 

.32729 

— 

.00013 

1.1 

- 

.19344 

- 

.21797 

.15225 

3.6 

— 

.32730 

— 

.00007 

1.2 

- 

.21442 

- 

.20098 

.18548 

3.7 

— 

.32731 

— 

.00004 

1.3 

- 

.23355 

- 

.18129 

.20629 

3.8 

- 

.32731 

— 

.0000? 

1.4 

- 

.25063 

- 

.1  6011 

.21550 

3.9 

— 

.32701 

— 

.00001 

1  .5 

- 

.26556 

- 

.13853 

.21 454 

4.0 

— 

.32731 

— 

.00000 

1.6 

- 

.27835 

- 

.11 747 

.20524 

4.1 

- 

.32731 

— 

.00000 

1.7 

- 

,?poi  n 

- 

.09768 

.18969 

4.? 

— 

.32731 

— 

.nnooo 

1.8 

- 

.29795 

- 

.07068 

.16995 

4.3 

— 

.32731 

— 

.O'-'nno 

1.9 

- 

.30510 

- 

.06377 

.14797 

4.4 

— 

.32701 

— 

.000^0 

2.0 

- 

.31078 

- 

.05010 

.12541 

4.3 

— 

.32701 

— 

.000*0 

2.1 

- 

■  .31520 

- 

.03866 

.10360 

4?6 

— 

.02701 

— 

.non** 

2.2 

- 

,31858 

- 

.02932 

.08349 

4.7 

— 

.02731 

— 

.ooo  np 

2.3 

- 

.32113 

- 

.0?189 

.06567 

4.8 

— 

.02731 

— 

.0*000 

2*4 

- 

,32301 

- 

*01610 

.04959 

4.9 

- 

.32731 

*00000 

-  .01758 

-  .01180 

-  .00778 
,On5o3 
.00320 
.00200 
.00122 
.00074 

-  .00043 

-  .00025 

-  .00014 

-  .00008 

-  .00004 

-  .00002 

-  .ooooi 

-  .oonoo 

-  .00000 
-  ,00000 
-  .00000 

-  .O00O* 

-  ,00000 

-  • 000PO 

-  .00000 

-  .ooooo 

-  *00000 


.04116 
.03012 
.02169 
.01515 
,01041 
,00700 
.00461 
.00298 
.00188 
.00116 
.00071 
.00042 
.00024 
.0001  4 
.00007 
.00004 
.00002 
.ooooi 

.opopo 
,  OOOOO 
.OOOOO 
.ooooo 
.  ooo  oo 

.OOOOO 

*00000 
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Equation  12 


X 

k 

k* 

kl 

* 

t 

.0 

.00000 

.onnno 

.58678 

2.5 

.3-2  70  7 

.00757 

*02835 

.1 

.00291 

.05798 

.56635 

2.6 

.32770 

.00515 

_ 

.02029 

.2 

.01146 

.11214 

.51209 

2.7 

.32812 

.00344 

— 

.01421 

.3 

.0251? 

.15960 

.43383 

2.8 

.32840 

.00226 

_ 

.00976 

.4 

.04310 

.19841 

.34050 

2.9 

.32859 

.00145 

— 

.00655 

.5 

.06447 

.22747 

.2401? 

3.o 

.32870 

.00091 

_ 

.00431 

.6 

.08825 

.24643 

.13966 

3.1 

.32878 

.00057 

_ 

.00278 

.7 

.11343 

.25559 

.04496 

3.2 

.3288? 

.00034 

_ ' 

.00175 

.8 

.13907 

.25577 

- 

.03941 

3.3 

.32885 

.00020 

— ► 

.00108 

.9 

.16433 

.24816 

- 

.11024 

3.4 

.32887 

.00012 

_ 

.00066 

l.n 

.18850 

.23424 

- 

.16561 

3.5 

.32888 

.00006 

.00039 

1.1 

.21102 

.21558 

- 

.20491 

3.6 

..32888 

.00003 

_ 

.00022 

1.2 

.23151 

.19377 

- 

.22867 

3.7 

.32888 

.00002 

— 

.00013 

1.3 

.24972 

.17031 

- 

.23833 

3.8 

.32889 

.00001 

.00007 

1.4 

.26556 

.14650 

- 

.23601 

3.9 

.32889 

.00000 

— 

.00004 

1.5 

.2  7905 

.12342 

- 

.22423 

4.0 

.32889 

.00000 

.00002 

1.6 

,  29°3r* 

.10189 

- 

.20560 

4.1 

.32889 

.  nooon 

.ononl 

1.7 

.29950 

.08244 

- 

.18269 

4.2 

.32889 

.  nfl«on 

.noono 

1.8 

.30687 

.06541 

- 

.15772 

4.3 

.32889 

.00000 

.00nno 

1.9 

.31266 

.05090 

- 

.13258 

4.4 

.32889 

.00000 

.noooo 

2.0 

.31713 

.03886 

- 

.10866 

4.5 

.32889 

.00000 

— 

.00000 

2.1 

.32051 

.02910 

- 

.08692 

4.6 

.32889 

.00000 

_ 

.00000 

2.2 

.32302 

.02138 

- 

.06793 

4.7 

.32889 

-  .00000 

— 

.0000O 

2.3 

.32485 

.01541 

- 

.05190 

4.8 

.32889 

-  .00000 

— 

,00000 

2.4 

,32615 

.01090 

- 

*03878 

4,9 

.32889 

-  ,00000 

- 

*00000 

Equation  13 

X  kfc 

5 

,0  .00000 

.1  -  .00057 

k5 

.00000 

.01148 

-  .11629 

-  .11219 

2.5 

2.6 

— 

.06250 

.06260 

— 

.00122 

.0008? 

.00477 

.00336 

.2 

— 

.00227 

— 

.02221 

-  .1012?  . 

2.7 

- 

.06267 

— 

.00054 

.00232 

.3 

— 

.00497 

- 

.03157 

-  .08531 

2.8 

- 

.06271 

— 

.00035 

.00156 

.4 

— 

.00852 

- 

.0391 7 

-  .06630 

2.9 

- 

.06274 

- 

.00022 

.00103 

.5 

— 

.01274 

- 

.04478 

-  .04587 

3.0 

— 

.06?  76 

— 

.OO013 

.00067 

.6 

— 

.01741 

- 

.04834 

-  .02550 

3.1 

— 

.06277 

— 

.00008 

.00042 

.7 

— 

.02234 

- 

.0499? 

-  .00644 

3.2 

— 

.06278 

— 

^nooo5 

.00026 

.8 

— 

.02734 

— 

.04971 

.01036 

3.3 

- 

.06278 

— 

.00002 

.00016 

.9 

— 

.03223 

- 

.04795 

.02425 

3.4 

— 

.06278 

— 

.0O001 

.0O009 

1.0 

— 

.03689 

— 

.04496 

.03488 

3.5 

- 

.06278 

—  . 

.00000 

.00005 

1.1 

— 

.04120 

- 

.04108 

.04216 

3.6 

- 

.06278 

— 

.00000 

.00003 

1.2 

— 

.04509 

- 

.03664 

.04626 

3.7 

— 

.06279 

— 

,00000 

.0000] 

1.3 

- 

.04852 

- 

.03193 

.0475?  • 

3.8 

— 

.06279 

— 

. 0  00  on 

,00001 

1  .4 

- 

.05147 

- 

.02721 

.04643 

3.9 

— 

.06279 

— 

.000OO 

.00000 

1  .5 

- 

.0  5397 

- 

.02270 

.04355 

4.0 

— 

.06279 

mm 

.OOOi"*0 

.00000 

1.6 

- 

.0560? 

- 

.01854 

.03941 

4.1 

— 

.06279 

— 

.noooo 

,00000 

1.7 

- 

.05769 

- 

.01484 

.03456 

4.2 

— 

.06279 

.OOOOO 

.OOOOO 

1.8 

- 

.0  5901 

- 

.01164 

.02944 

4.3 

— 

.06279 

— 

.00000 

,00000 

1.9 

- 

.O6OO4 

- 

.00805 

.02441 

4.4 

— 

.06279 

.onoon 

,00000 

2.0 

- 

.0608? 

- 

.00674 

.Ol 97? 

4.5 

— 

.06?  79 

.00000 

.00000 

2.1 

- 

.06140 

- 

.00498 ■ 

.01 555 

4.6 

— 

•w06?79 

, 00000 

.OOOOO 

2.2 

- 

.061  83 

- 

.00361 

.01197 

4.7 

— 

.06279 

.00000 

.OOOOO 

2.3 

- 

.06213 

- 

.00257 

.00901 

4.8 

— 

.06279 

.00000 

.OOOOO 

2.4 

— 

.06235 

— 

.00179 

.00663 

4,9 

- 

.06279 

,00000 

*  OOOOO 

Equation  lV- 

X, 

h  , 

h* 

hV 

* 

V 

.0 

.OOOOO 

.oonoo 

. 

.72468 

.1 

. 00360 

.07190 

.70791 

.2 

.01427 

.14059 

.66139 

.3 

.03153 

.20337 

.59078 

.4 

,05468 

.25813 

.50176 

.5 

.08283 

.30330 

.39995 

.6 

,11498 

.33788 

.29088 

.7 

,15004 

.36141 

.17982 

.8 

.1 8690 

.  37394 

.07172 

.9 

.22448 

.37600 

— 

.02899 

1.0 

.26178 

.36851 

— 

.11858 

1.1 

.29791 

.35275 

- 

.19409 

1.2 

.33211 

.33023 

— 

.25354 

1.3 

.36379 

.30261 

- 

.29597 

1.4 

. 39252 

.271-60 

— 

.32145 

1.5 

.41805 

.23885 

— 

.33104 

1.6 

.44028 

.20586 

- 

.32660 

1.7 

.45926 

.17392 

— 

.31059 

1.8 

.47514 

.14403 

— 

.28584 

1.9 

.48816 

.11694 

— 

.25527 

2.0 

.49864 

.09308 

- 

.22162 

2.1 

.50689 

.07264 

— 

.18733 

2.2 

,51328 

.05557 

— 

.15432 

2.3 

.51811 

.04168 

— 

.12400 

2.4 

.52171 

,03065 

— 

.09726 

Equation  15 

h" 

\ 

h 

h' 

5 

5 

5 

.0  ■ 

.00000 

.oonoo 

.50433 

.1 

.00250 

. °4973 

.48391 

.2 

.0098? 

.09572 

.43168 

.3 

.02143 

.13541 

.35961 

.4 

.03664 

.16731 

.27750 

.5 

.05462 

.19083 

.19306 

.6 

.0  7453 

.20604 

.11202 

.7 

.09557 

.21349 

.03838 

.8 

. 11700 

.21405 

- 

.02535 

.9 

,13818 

.20879 

- 

.07792 

1  .0 

.1 5860 

.  19885 

- 

.11900 

1.1 

. 17784 

.18536 

— 

.14895 

1.2 

.1 9559 

.16940 

- 

.16857 

1.3 

.21167 

.15195 

- 

.1  7899 

1.4 

.22596 

.1  3387 

- 

.18147 

1.5 

.23845 

.11587 

- 

.17738 

1.6 

.24916 

.09856 

- 

.16809 

1.7 

.25820 

.08238 

- 

.15494 

1.8 

.26569 

.06766 

- 

.13921 

1.9 

.27178 

.05459 

- 

.12208 

n  o 

<_  • 

.27666 

.04325 

- 

.10457 

2.1 

.28049 

.03366 

- 

.  0  8  7  61 4 

2.2 

.28345 

.n?C7l 

- 

.071 64 

2.3 

.28569 

.01  9?.7 

- 

.05733 

2,4 

, 2  87  3  5 

,ol 41 8 

- 

,04486 

2,5- 

,52433 

.02210 

— 

,07448 

2.6 

,52620 

.01562 

— 

.06574 

2.7 

.52751 

.01083 

- 

.04077 

2.8 

.52841 

.00736 

- 

.02916 

2.9 

.52902 

.00490 

- 

.02040 

3.0 

.52942 

.00320 

- 

.01396 

3.1 

.52967 

.00205 

- 

.00934 

3.2 

.52984 

.00128 

— 

.00612 

3.3 

.52994 

.00079 

— 

.00393 

3.4 

.53000 

.00047 

- 

.00246 

3.5 

.53004 

.00028 

— 

.00151 

3.6 

. .53006 

.00016 

— 

,00091 

3.7 

.53008 

.00009 

— 

.00053 

3.8 

.53008 

.00005 

- 

.00031 

3.9 

.53009 

.00002 

— 

.OOOl 7 

4.0 

.53009 

•ooooi 

— 

.00009 

4.1 

.53009 

.oonoo 

— 

.00005 

4.2 

.53009 

.oonoo 

— 

.00002 

4.3 

.53009 

.00000 

— 

.OOOOI 

4.4 

.53009 

.00000 

— 

.OOOOO 

4.5 

.53009 

.00000 

- 

.000OO 

4.6 

.53009 

.00000 

— 

.ooooo 

4.7 

.53009 

.ooooo 

- 

.OOOOO 

4.8 

.53009 

, oonoo 

- 

.ooooo 

4,9 

.53009 

.oonoo 

- 

*00000 

2.5 

.28856 

.01024 

— 

*03434 

2.6 

.28943 

.00725 

— 

.02571 

2.7 

.29004 

.00503 

— 

.01 884 

2.8 

.29046 

.00343 

— 

.01 350 

2.9 

.29074 

.00220 

— 

.00947 

3.0 

.29093 

.00150 

- 

.00650 

3.1 

.29105 

.00096 

— 

.00438 

3.2 

.29113 

.00060 

— 

.00287 

3.3 

.29118 

.00037 

— 

.00185 

3.4 

.29121 

.00022 

— 

.00116 

3.5 

.29122 

.00013 

— 

.00071 

3.6 

.291 23 

.00007 

- 

.00043 

3.7 

.291.24 

.00004 

— 

.00025 

3.8 

.291 24 

,00002 

— 

.00014 

3.9 

.29125 

. ooooi 

— 

, nnnnp 

4.0 

.29125 

.Ooo^o 

— 

.00004 

4.1 

.29125 

.ooooo 

— 

,00002 

4.2 

.29125 

.00000 

— 

,nonol 

4.3 

.29125 

.0O0OO 

— 

.OOnon 

4.4 

.29125 

.OOOOO 

- 

.ooooo 

4.6 

.29126 

.ooooo 

— 

, nonon 

4.6 

.29126 

.OOOOO 

— 

.  nonon 

4.7 

.29125 

.  oonoo 

— 

, oonoo 

4.8 

.29125 

.OOOOO 

— 

,  noooo 

4,o 

.201 26 

-  .ooooo 

— 

.nonon 
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Equation  1 6 


\ 

h6 

h6 

h6 

.0 

#oonon 

. onnoo 

— 

.56841 

.1 

- 

.00283’ 

- 

.05641 

- 

.55580 

•  2 

- 

.01120 

- 

.11040 

— 

.52067 

.3 

- 

.02476 

— 

.15993 

— 

.46711 

.4 

- 

.04298 

- 

.20335 

- 

.39921 

.5 

- 

.06519 

— 

.23943 

- 

.32107 

.6 

- 

.09060 

- 

.26736 

- 

.23678 

.7 

- 

.11837 

- 

.28672 

- 

.15032 

.8 

- 

.14765 

- 

.29748 

— 

.06546 

.9 

- 

.  17759 

- 

.29998 

.01432 

1.0 

- 

.20740 

— 

.29488 

.08600 

1.1 

- 

.23635 

- 

.2831 3 

.14713 

1.2 

- 

.26384 

— 

.26587 

.19598 

1.3 

— 

.2  8938 

- 

.24438 

.231 57 

1.4 

- 

.31262 

— 

.22000 

.25378 

1.5 

- 

.33333 

— 

.19404 

.26325 

1.6 

- 

.35142 

- 

.  16773 

.26131 

1.7 

- 

.36690 

- 

.14210 

.24984 

1.8 

- 

.37989 

— 

.11800 

.23104 

1.9 

- 

.39057 

- 

♦09605 

.20723 

2.0 

- 

.39919 

— 

.07665 

.18063 

2.1 

- 

.40600 

- 

.05995 

.15323 

2.2 

- 

.41127 

- 

.04597 

.12664 

2.3 

- 

.41528 

— 

.03456 

.10207 

2.4 

- 

.41826 

- 

,02547 

.08028 

Equation  17 


\ 

hT 

.0 

. OOOO" 

,  OO  HAf) 

,26501 

.1 

. 00132 

.02637 

.26122 

.2 

.00525 

."5199 

.25006 

.3 

.01167 

.07615 

.23191 

•  4 

.02040 

.09816 

.20733 

.5 

.03121 

.11743 

.17710 

.6 

.04378 

.  13343 

.14228 

.7 

.05778 

.14577 

.10416 

.8 

.0  72  81 

.15420 

.06424 

.9 

.0  8848 

.15862 

.02418 

1.0 

.1^440  . 

. 1 5909 

- 

.01431 

1.1 

.12018 

.1558  6 

- 

.04962 

1.2 

. 1 3546 

.14932 

- 

.08029 

1.3 

. 14995 

.  1  4000 

- 

.10518 

1  .4 

. 1 6339 

. 1 285" 

- 

♦  1  2  3  5  ^ 

1.5 

. 1 7560 

.  H  551 

— 

.13500 

1.6 

. 1 8647 

.10170 

- 

.13999 

1.7 

.19594 

.08772 

- 

.13884 

1.8 

.  2  °4"  3 

.0741  1 

- 

.13257 

1.9 

.21079 

.06133 

- 

.12236 

2  •" 

,21633 

.04972 

- 

.10945 

2.1 

.22078 

.03949 

- 

#no50A 

2.2 

.22428 

.0307? 

- 

.08028 

2.3' 

.22698 

.02742 

- 

.06600 

2.4 

.22901 

.01749 

- 

.0^287 

2.5  . 

- 

.42044 

.  — 

,01840 

*06165 

2.6 

- 

.42200 

- 

.01303 

.04625 

2.7 

- 

.42309 

- 

.00904 

.03391 

2.8 

- 

.42384 

- 

.00616 

.02430 

2.9 

- 

.42435 

- 

.00411 

.01703 

3.0 

- 

.42469 

- 

.00269 

•  .01168 

3.1 

- 

.42490 

- 

.00172 

.00783 

3.2 

- 

.42504 

- 

.00108 

.00514 

3.3 

- 

.42513 

- 

.00066 

.00330 

3.4 

- 

.42518 

- 

.00040 

.00707 

3.5 

— 

.42521 

- 

.00023 

.00127 

3.6 

- 

.42523 

- 

.00013 

,00077 

3.7 

- 

..42524 

- 

.00007 

,00045 

3.8 

- 

.42525 

- 

.00004 

.00026 

3.9 

- 

.42525 

- 

.00002 

.00014 

4.0 

- 

.42525 

- 

.00001 

.00008 

4.1 

- 

.42525 

- 

.OOOOO 

.00004 

4.2 

- 

.42525 

- 

.00000 

.00002 

4.3 

- 

.42525 

- 

.00000 

.00001 

4.4 

- 

.42525 

— 

.00000 

.00000 

4.5 

- 

.42525 

- 

.00000 

.00000 

4.6 

- 

.42525 

- 

.00000 

.00000 

4.7 

- 

.42525 

- 

.00000 

,00000 

4.8 

- 

.42525 

— 

.00000 

.00000 

4.9 

- 

.42525 

- 

.00000 

,00000 

2.5 

,23051 

.01279 

— 

.041 3o 

2.6 

.23160 

.00917 

— 

.03148 

2.7 

.23238 

.00643 

— 

.02342 

2.8 

.23292 

.0044? 

- 

.01 70? 

2.9 

.23328 

-.002  98 

— 

.01208 

3.0 

.23353 

.00197 

— 

.00838 

3.1 

.23369 

.00127 

— 

.00568 

3.2 

.23379 

.00081 

— 

.00377 

3.3 

.23386 

.OOO50 

— 

.00744 

3.4 

.23390 

.00030 

— 

.00155 

3.5 

.23392 

.00018 

— 

.00096 

3.6 

.23393 

,00010 

— 

.  0On  5  p 

3.7 

.23394 

.00006 

— 

.OOO34 

3.8 

.23395 

.00003 

- 

,nnn?n 

3.9 

.23395 

.OOOOj 

— 

.00011 

4.0 

.23395 

.00001 

— 

.ooon^ 

4.1 

.23395 

.00000 

- 

.00003 

4.2 

.23395 

.0O00O 

— 

.oooni 

4.3 

.23395 

,00000 

#nnnnn 

4.4 

.23395 

.  oonnn 

— 

.onono 

4.e; 

.23395 

.noooo 

— 

. OOOOn 

4.6 

.237Q5 

.00000 

— 

. oonnn 

4.7 

.23395 

.000^0 

— 

. onnnn 

4.8 

.23395 

.onnnn 

— 

. ooonn 

4.9 

.2339^ 

#nnooo 

• 

.onnnn 
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Equation  18 

X  h8  h8  h8 


•  0 

,00000 

.ooooo 

- 

.51085 

2.5 

— 

.‘40874 

— 

.02036 

,06682 

.1 

- 

.00254 

- 

.05076 

- 

.50128 

2.6 

— 

.41 047 

— 

.01452 

.050^8 

.2 

- 

.01009 

- 

.09967 

- 

.474-15 

2.7 

- 

.41169 

- 

.01015 

.03739 

.3 

- 

.02236 

- 

.14509 

- 

.43191 

2.8 

— 

.41254 

- 

.00695 

.02701 

V4  . 

*.03894 

- 

.18563 

- 

.37711 

2.9 

- 

.41311 

— 

.00466 

.01 907 

.5 

- 

.  03-92  9 

- 

.22018 

— 

.31249 

3.0 

- 

.41349 

- 

.00307 

.01316 

.6 

— 

.08275 

- 

.24790 

- 

.24093 

3.1 

— 

.41 374 

- 

.00198 

.00889 

.7 

- 

.1 0862 

- 

.26824 

- 

.16545 

3.2 

— 

.41390 

— 

.00125 

.00587 

.8 

— 

.13618 

- 

.28096 

- 

.08918 

3.3 

- 

.41  400 

- 

.00077 

.00379 

.9 

— 

.16456 

.2861 5 

- 

.01521 

3.4 

— 

.41406 

— 

.00047 

.00289 

1.0 

- 

.19314 

- 

.2841 8 

.05354 

3.5 

— 

.41410 

— 

.00028 

.00148 

1.1 

- 

.22118 

- 

.27570 

.11446 

3.6 

- 

.41412 

- 

.00016 

.00089 

1 .2 

- 

. 24809 

— 

.26162 

.16542 

3.7 

— 

.41413 

— 

.000^9 

.00053 

1.3 

— 

.27336 

- 

.24300 

.20491 

3.8 

— 

.41414 

- 

.00005 

.00030 

1.4 

- 

.29658 

- 

.22104 

.23217 

3.9 

— 

.41414 

— 

.00002 

.00017 

1.5 

- 

,31750 

- 

.19697 

.24724 

4.0 

— 

.41415 

— 

.00001 

.00009 

1.6 

- 

.33595 

- 

.17198 

.25087 

4.1 

— 

.41415 

— 

.00000 

.00005 

1.7 

- 

.35190 

- 

.1471  3 

.24449 

4.2 

- 

.41415 

- 

.00000 

,00002 

1.8 

- 

.36541 

- 

.12335 

.22997 

4.3 

— 

.41415 

— 

.0O000 

.O00O1 

1.9 

- 

.3  7663 

- 

.10134 

.20946 

4.4 

— 

.41415 

— 

.00000 

.0^000 

2.0 

- 

.38576 

- 

.081 58 

.18515 

4.5 

— 

.41415 

— 

.00000 

,  0000  r' 

2.1 

- 

.39303 

- 

.06436 

.1  5911 

4.6 

- 

.41415 

— 

.00000 

.ooooo 

2.2 

- 

. 39872 

— 

.04976 

.13308 

4.7 

— 

.41415 

— 

.ooooo 

.ooooo 

2.3 

— 

.40307 

- 

.03770 

.10846 

4.8 

— 

.4141 5 

— 

.00000 

.ooooo 

2.4 

- 

,40634 

- 

.02799 

.08618 

4.9 

- 

.41415 

- 

,00000 

40oono 

Equation  19 


X 

h 

hi 

h" 

9 

9 

9 

.0O00O 

.00^00 

- 

.95816 

2.5 

- 

.  8  5.0  7  7 

— 

,04839 

.15585 

.1 

- 

. O0478 

- 

.09537 

- 

.94503 

2.6 

— 

.85440 

— 

.03471 

.11 8°7 

.2 

- 

.01898 

- 

.1881 5 

- 

.90631 

2.7 

— 

.85733 

— 

.02439 

.08857 

.3 

- 

.04224 

- 

.27582 

- 

.84322 

2.8 

— 

.85.937 

— 

.01679 

.06441 

.4 

- 

.07390 

- 

.35604 

- 

.75751 

2.9 

- 

.86076 

— 

.01132 

.04576 

.5 

- 

,1131 2 

- 

.42665 

- 

.65159 

3.0 

- 

.86169 

— 

.00748 

.03177 

.6 

- 

.1 5885 

- 

.48579 

- 

.52876 

3.1 

- 

.86729 

— 

.00485 

.921  ^ 

.7  • 

- 

.20985 

- 

. 53198 

- 

,3937° 

3.2 

- 

.86268 

— 

.00807 

.01431 

.8 

- 

.26478 

- 

.5641 8 

- 

.24996 

3.3 

- 

.86293 

- 

.00191 

•  009?  9 

.9 

- 

.32221 

- 

.58190 

- 

.10477 

3.4 

- 

.86308 

— 

.00116 

.00589 

1.0 

- 

.38068 

- 

.6fl527 

.03626 

3.5 

— 

.86317 

— 

.00069 

.00366 

1.1 

- 

.43881 

- 

.57499 

.16704 

3.6 

- 

.86323 

— 

,00040 

.00227 

1.2 

- 

.49527 

- 

. 55239 

.28201 

3.7 

— 

,86326 

— 

,ooo73 

.non? 

1.3 

- 

.54893 

- 

.61 027 

.37661 

3.8 

— 

.86327 

— 

.00013 

.00077 

1.4 

- 

.598^5 

- 

.47785 

.44766 

3.9 

— 

.86328 

— 

,00007 

.00044 

1.5 

- 

.64431 

- 

.43057 

.4936^ 

4.0 

- 

.86329 

— 

.OOO03 

,ooo?4 

1  .6 

- 

.68486 

- 

.37995 

.51490 

4.1 

- 

.86329 

— 

.00002 

,00013 

1.7 

- 

.72027 

- 

.32835 

.51335 

4.7 

- 

.86329 

— 

.OOOOl  ' 

.00007 

1.8 

- 

.750  57 

- 

.27793 

.49230 

4.3 

— 

.86330 

— 

.00000 

.00003 

1.9 

- 

.7759^ 

- 

.73040 

.46600 

4.4 

— 

.86330 

— 

.00000 

.  oono  1 

2.0 

- 

.79679 

- 

.1  8708 

.40910 

4.5 

— 

.86330 

— 

. ooooo 

. 0 00  00 

2.1 

- 

.81354 

- 

.14878 

.36621 

4.6 

— 

.86330 

— 

.OOnno 

.  0  n  n  r>  r\ 

2.2 

- 

.82673 

- 

.  11  600 

.30149 

4.7 

— 

.86330 

— 

.OOOOO 

.  A  o  r>  r\  n 

2,3 

- 

.83690 

- 

.08843 

.24833 

4.8 

— 

.86330 

— 

.00000 

O  r>  0  0  r\ 

. 

2.4 

- 

. 84458 

- 

.06609 

.19976 

4.0 

- 

,R633o 

- 

.00000 

,  nooon 

57 


